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Abstract 



The scale invariance of the quantum fluctuations in de Sitter space leads to the appearance 
of de Sitter symmetry breaking infra-red logarithms in the graviton propagator. We investi- 
gate physical effects of soft gravitons on the local dynamics of matter fields well inside the 
;2j ! cosmological horizon. We show that the IR logarithms do not spoil Lorentz invariance in 

^.' scalar and Dirac field theory. The leading IR logarithms can be absorbed by a time depen- 

dent wave function renormalization factor in the both cases. In the interacting field theory 
with A(^^ and Yukawa interaction, we find that the couplings become time dependent with 
definite scaling exponents. We argue that the relative scaling exponents of the couplings are 
gauge invariant and physical as we can use the evolution of a coupling as a physical time. 
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1 Introduction 



In de Sitter (dS) space, the degrees of freedom outside the cosmological horizon increase 
with cosmic evolution. This increase leads to the de Sitter symmetry breaking term in the 
propagator of a massless and minimally coupled scalar field. It is a direct consequence of the 
scale invariant fluctuation spectrum [11 [2l [3j. So in some field theoretic models in dS space, 
physical quantities may become time dependent through the propagator. The gravitational 
field on dS background is a candidate which induces such infra-red (IR) effects. It is because 
the gravitational field is massless and contains minimally coupled modes [4J . 

In order to investigate interacting field theories in dS space, we need to employ the Schwinger- 
Keldysh perturbation theory [5l [6] . The IR effects at each order manifest as the polynomials 
in the logarithm of the scale factor of the universe: log a(r) \J\. In an interacting scalar 
field theory with polynomial interactions such as Ay?^ theory [S], this IR logarithms in the 
propagators give rise to powers of IR logarithms in the amplitudes. Let us consider the 
expectation value of the energy-momentum tensor for example. The leading power of IR 
logarithms is given by the number of the propagators of the diagram. They make an effective 
cosmological constant time dependent and thus break de Sitter symmetry. Therefore we need 
to sum these leading IR logarithms to understand the long term evolution as the effect of 
IR logarithms becomes large if we wait long enough. 

Remarkably a simple physical picture holds in the leading IR log approximation as follows. 
A scalar field is not completely frozen beyond the horizon scale as it is constantly jolted by 
the modes coming out of the horizon. A scalar field performs a random walk in the scalar 
field space (1 dimension) which is consistent with the linear growth of the propagator with 
respect to cosmic time t at the initial stage. Eventually it reaches an equilibrium state in a 
potential well as de Sitter symmetry is restored [Qj, |10j . We have investigated IR logarithms 
in the non-linear sigma models which contain massless minimally coupled scalar fields with 
derivative interactions. In the expectation value of the energy-momentum tensor, we have 
shown that the leading IR logarithms cancel to all orders [TT] . 

It is important to understand the IR effects in quantum gravity as the propagator of gravitons 
contains IR logarithms in de Sitter space. The case is very strong here as gravitons exist 
in our Universe which is of de Sitter type. In this paper we investigate physical effects 
of soft gravitons on microscopic physics in de Sitter space. We focus on the dynamics of 
matter fields with sub-horizon momentum scale as it is directly observable. In contrast, the 
momentum scale of matter fields is taken to be super- horizon scale in ^^ [T3] . We find 
that the super-horizon gravitons influence the matter field dynamics inside the cosmological 
horizon. Our predictions are thus verifiable by direct observations in principle. 

In investigating the soft graviton effects to local matter field dynamics, it is a non-trivial 
question whether they preserve the Lorentz invariance. To verify the Lorentz invariance, we 
investigate the kinetic terms of matter fields, specifically the relative weight between the time 
derivative term and the spatial derivative term. In interacting field theories we investigate 
how the IR effects influence their coupling constants. As specific examples, we adopt Xcf)^ 
and Yukawa theory with dimensionless couplings . 



^We recall here the fractal dhxiension of random walk is 2. 



The organization of this paper is as follows. In Section 2, we quantize gravitational field on 
dS background. We identify the graviton modes which exhibit IR logarithm. In Section 3, 
we evaluate the quantum equation of motion with respect to matter fields which are dressed 
by soft gravitons at the one- loop level. First, we adopt free field theories and investigate 
whether the Lorentz invariance is preserved. Second, we evaluate the effective coupling 
constants in 0^ and Yukawa theory. We find that the IR effects from gravitons preserve 
the Lorentz invariance. The effective coupling constants are found to decrease with cosmic 
expansion with definite scaling exponents. In Section 4, we vary the gauge parameter of the 
graviton propagator to investigate the gauge dependences of the results obtained in Section 
3. We show that the relative scaling exponents of the couplings are gauge invariant and 
observable. We conclude with discussions in Section 5. 



2 Gravitational field in de Sitter space 



In this section, we compute the gravitation propagator in de Sitter (dS) space. In the 
Poincare coordinate, the metric in dS space is 

ds^ = -dt^ + a^{t)dx^, a{t) = e^\ (2.1) 

where the dimension of dS space is taken as -D = 4 and H is the Hubble constant. In the 
conformally flat coordinate, 

{gi,u)AS = a'^{r)r]f,^, a{T) = - — . (2.2) 

Here the conformal time r (— oo < r < 0) is related to the cosmic time t as r = — -^e^^*. 
We assume that de Sitter space begins at an initial time tj with a finite spacial extension. 
After a sufficient exponential expansion, the de Sitter space is well described locally by the 
above metric irrespective of the spacial topology. The metric is invariant under the scaling 
transformation. 

r -» CT, x" -)> cx\ (2.3) 

It is a part of the S'0( 1,4) de Sitter symmetry. 

In dealing with the quantum fluctuation whose background is dS space, we adopt the fol- 
lowing parametrization: 

g^, = n\xyg,,, n{x) = a{T)e^^^-\ (2.4) 

det^^, = -l, ^^. = (e'''^(^))^,, (2.5) 

where k is defined by the Newton's constant G as k^ = IGvrG. To satisfy (12. 5p . h^y is traceless 

v'^'^Ku = 0. (2.6) 



By using this parametrization, the components of the Einstein action are written as follows. 
We keep a parameter D to specify the dimension for generality: 

v/=^ = f]^, (2.7) 

R = n-^R - 2{D - l)n-^g^"'Vf,d^n -{D-1){D- A)n-^g'"'d^ndM, (2.8) 

where R is the Ricci scalar constructed from g^^iy 

R = -d,d,r - \rrr^d,~9p.d,~9.fi + ^g^'^rr^d^g.^Qp-guf^- (2.9) 

From (12. 7p and (12. 8p , the lagrangian of gravity is 

gravity = ^V^[i? - iP - l)(i^ " 2)1/2] ^2.10) 

= 4 l^^^^R + (D-i)(D- 2)n^~^~g^'^dsid^n -{d-imd- 2)H^n^] . 

Note that the lagrangian density is defined including y/—g in this paper. 

In order to fix the gauge degrees from general coordinate invariance 

x'^ = Xf, + e^, (2.11) 

dfMu 9^J.l^ QfipOu^ QupOjiS Opg^i^E , 

we introduce the gauge fixing term [1] 

/:gf = —a^'-VF^F,, (2.12) 

Ff^ = dph; -{D- 2)d^w + {D- 2)/i/9,loga + 2{D - 2)wd^\oga. 

The corresponding ghost term at the quadratic level is 

^Cghost = - a'''^dj>^r,''''{r,ppd, + r^^pd^ + 2r,p,dpi\oga)]h'' (2.13) 

where h^' is a ghost field and ¥ is an anti-ghost field. From (EHD, ( 1^^ - f ETT^ and (EUSJ, 
the quadratic part of the total Lagrangian density is 

-^quadratic = a'"'^ {]-D{D - 2)r,>^'' d,wd,w - ^.^ ^y dph'>'>d,h''> - ]r]>^''d,h' ^d^h^ 
Z A[1J — Ij 4 

+ a^i/2{ - 2{D - 2)w' + 2{D - 2)w/i°° - \d - 2)/i°°/i°° (2.14) 

+ -(D - 2)/i0*/i0* + (Z^ - 2)6°60}. 



Here we have decomposed /i* •, i,j = l,---,D — l into a trace part and a traceless part 

^^• = ^^• + D^^h\s^^ = h^^ + ^/^°°5v (2.15) 

fl2.14p is diagonalized as 

-^quadratic = a"" [ ^a-'^v'^'^d^Xd^X - -^a-^i^^'^ dji' ^d,h\ - a-^T^^^ d,Vd,V (2.16) 

+ -a-^r]^''d^h'^'d^h'^' + -{D- 2)H^h°'h^' 



1 

a 

2 



^r]''''d^Yd^Y -{D- 3)H^Y^ 



where X, Y are 



X = {D- 2)J^-^w - ^ =/i°°, Y = J /^ ^ (/i°° - 2«;). (2.17) 

^ ^V/^-3 ^p_i)(D_3) y2(D-3)^ ; I ; 

f l2.16p contains three types of fields, massless and minimally coupled fields: X, /i* •, 6\ 6*, 
massless conformally coupled fields: /i°*,6°,6°, and Y. 

We restrict to the D = A case in the subsequent discussion. In this case, Y corresponds with 
a massless conformally coupled field. We list gravitation propagators as below 



{X{x)X{x')) = -{^{xMx')), (2.18) 

2 
— ( 
3 



Ch'^{x)h',{x')) = i6%i + 6\6^' - ^6^^^S\){^{xMx')) 



{b\x)V {x')) = 6'^ip{xMx')), 

(/i°'(x)/i°^(x')) = -<5'^(0(s)0(x')), (2.19) 

{Y{x)Y{x')) = (0(x)0(x')), 
{b%x)W{x')) = -(0(x)0(x')), 

Here ip denotes a massless and minimally coupled scalar field and denotes a massless 
conformally coupled scalar field 

{^{xMx')) = §^{- - ^ logy + ^ log a(r)a(r') + 1 -7}, (2.20) 

(0(x)0(x')) = f-,-, (2.21) 

where 7 is Euler's constant and y is the dS invariant distance 

y = . (2.22) 



It should be noted that the propagator for a massless and minimally coupled scalar field has 
the dS symmetry breaking logarithmic term: log a(r)a(r'). To explain what causes the dS 
symmetry breaking, we recall the wave function for a massless and minimally coupled field 

0p(x) = ^(l-z^)e---^-- (2.23) 

Inside the cosmological horizon P = p/a^r) ^ iJ <^ p|r| ^ 1, this wave function approaches 
to that in Minkowski space up to a cosmic scale factor 

Hr 

0p(x) ~ -^e-^^+P-^ (2.24) 



On the other hand, the behavior outside the cosmological horizon P ^ iJ is 

0p(x) ~ ^e^^-. (2.25) 

The IR behavior indicates that the corresponding propagator has a scale invariant spectrum. 
As a direct consequence of it, the propagator has a logarithmic divergence from the IR 
contributions in the infinite volume limit. 

To regularize the IR divergence, we introduce an IR cut-off Eq which fixes the minimum value 
of the comoving momentum 

dP. (2.26) 

With this prescription, more degrees of freedom (d.o.f.) go out the cosmological horizon with 
cosmic evolution. Due to the increase, the propagator acquires the growing time dependence 
which spoils the dS symmetry. In tribute to its origin, we call the dS symmetry breaking 
term the IR logarithm. Physically speaking, I/eq is recognized as an initial size of universe 
when the exponential expanding starts. For simplicity, we set Eq = H in (12.201) . 

As there is explicit time dependence in the propagator, physical quantities can acquire time 
dependence through the quantum loop corrections. We call them the quantum IR effects in 
dS space. By the power counting of the IR logarithms in quantum gravity, the leading IR 
effects at ra-loop level is estimated as (k^ if ^ log a(r))". It indicates that even if k^H^ <^ 1, 
the quantum effects can eventually grow up to the tree level magnitude . This is the reason 
why we focus on the quantum IR effects in dS space. Before concluding this section, we 
introduce an approximation. Focusing on the IR effects, we can neglect conformally coupled 
modes of gravity since they do not induce the IR logarithm. By applying this approximation, 
we can identify the following two modes as 

^/3 
/i°o ~ 2u; ~ —X. (2.27) 



3 Quantum equation of motion 



In the preceding section, we have reviewed that gravitational field contains massless and 
minimally coupled modes and the corresponding propagator is time dependent due to the 



increase of d.o.f. outside the cosinological horizon. In this section, we investigate how the 
quantum IR effects from gravitons influence the local dynamics of the matter fields. More 
specifically, we evaluate the quantum equation of motion for matter fields including the 
quantum IR effects from gravitons. 

To begin with, we review how to derive the quantum equation of motion on a time dependent 
background p^. Let us represent the vacuum at t — )■ — cxd as |m), and t — )■ oo as \out). In 
the Feynman-Dyson formalism on a fiat background, it is presumed that \out) is equal to 
\in) up to a phase factor. On the other hand, we can't prefix \out) in de Sitter space. The 
correct strategy is to evaluate vacuum expectation values (vev) with respect to \in): 

{Oh{x)) = (m|Tc[f/(-oo,oo)f/(oo,-oo)C/(x)]|m). (3.1) 

where Oh and Oj denote the operators in the Heisenberg and the interaction pictures re- 
spectively. U{ti,t2) is the time translation operator in the interaction picture 

rh 

U{ti,t2) = exp {i / d^x 6Ci{x)}. (3.2) 

Jt2 

6£ denotes the interaction term of the lagrangian. It is crucial that the operator ordering 
Tc specified by the following path is adopted here 

c 

-oo ^ X > +CO , (3.3) 
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l-oo 


J-oo 





We call it the Schwinger-Keldysh formalism. Since there are two time indices +, — in this 
formalism, the propagator has four components 

(^+(x)vp+(x')) {v+{x)v-{x'))\^({Tv{x)v{x')) {^{x')^{x))\ 

{^4x)M^')) {^4x)^.{x'))J [ {^{xMx')) {f^{xMx'))J ' ^^-^^ 

where T denotes the anti-time ordering. 

We introduce the external source J+, J_ for each path and evaluate 

Z[J^, J_] = (m|Tc[f/(— oo, oo)?7(oo, — oo) exp {i / d'^x (J+(/9+ — J_(/9_)}]|m). (3.5) 

The generating functional for the connected Green's functions is 

zW[J+,J^] = log Z[J+,J^]. (3.6) 

We define the classical field as 

M^) = CAB ^Z^j^ V ^ A5 = +,-, (3.7) 

dJB{x) 



Cab 



1 
-1 



(3. 



By taking the limit J+ = J_ = J in (13. 7p . we obtain the vev of ip where the action contains 
the additional J(f term 



{^ix))\jv> = <^+(a;)|j+=j_=j = (^-(x)|j+=j_=j. 
Finally, we turn off the source term J = 



(^(x)) 



SW[J+, J_ 



6J-^-{x) 



5W[J+, J_ 



J+=J_=0 



6J4x) 



J+=J_=0 



(3.9) 



(3.10) 



The effective action is obtained after the Legendre transformation 



r[^+,^_] = w[j+, j_] - jd^'x (j+^+ - j_^_), 



where J+^_ are given by 0+,- as follows 

Ja(x) = 



-Cab 



5T[ip+,if^ 



5(pb{,x) 
From (13. 9 p (I3.12p . we obtain in the limit (^+ = (p^ = (p 



J{x) 



bYYp^,(p. 



#+(x) 



bYYp^,(p. 



^+=9-=^ 



6ip^{x) 



v+=f-='p 



(3.11) 



(3.12) 



(3.13) 



In the absence of the external source, the exact equation of motion is obtained including 
quantum effects 



5T[0+,0. 



60+{x) 



5T[ip+,0_ 



f+=v-=<f 



5ip-{x) 



0. 



(3.14) 



ip+=ip-=ip 



We call these identities quantum equation of motion in this paper. In the following subsec- 
tions, we evaluate these identities in concrete models to understand quantum IR effects from 
gravitons to matter field dynamics at the one-loop level. 



3.1 Free field theories 



Let us investigate the effects of soft gravitons on the local dynamics of a free scalar and Dirac 
field. First we investigate a massless conformally coupled scalar field. The corresponding 
action is 



S = I ./^d'x [ - \g^''d,<pd,<l) - ^R<p'] . 



(3.15) 



For convenience, we redefine the matter field Q(j) ~^ 0- 

S = j d'x[- hi^^d.^d,^ - ^^02]. (3.16) 

To obtain the quantum equation of the matter field, we decompose it into the classical field 
and the quantum fluctuation 

0->(^ + 0. (3.17) 

By differentiating 03.161) with respect to 0+ as in 03.141) . we obtain the quantum equation of 
motion. The quantum equation of motion up to the one-loop level is expressed as 

d^{x) + \H'd,{lih^^)^{x)(h;)^{x))d;^{x)) (3.18) 

where d'^ = rj^'^d^dv and we have taken the limit </>+ = 0_ = (j). Furthermore, we have 
neglected the contribution from differentiated gravitational field. It is because we focus on 
the IR logarithms and the differentiated gravitational field doesn't induce them. Since we 
consider the dynamics of matter fields whose momentum scale is larger than the Hubble 
scale, there is no loss of generality to assume that the scalar field is conformally coupled to 
gravity. 

From 02.18P and (I2.27p . we have only to focus on the following propagators to extract massless 
and minimally coupled modes from gravity 

(/iO°(x)/i°°(a;')) ^ --{^{x)^{x')), (3.19) 



1 



{h!'\x)h%{x')) ~ —JMx)^{x')) 



4" ^ 



3 



{W ^{x)h\{x')) ^ {6''6,i + 5\5.^ - -5^^5\){v{x)v{x')). 

By adopting an ultra-violet (UV) regularization, the propagator at coincident point is esti- 
mated as follows 

{ip{x)ip{x)) = (UV divergent const) + -— -loga(r). (3.20) 

47r^ 

Henceforth we focus on the time dependent dS symmetry breaking part 

{ifixMx)) ^ ^\oga{T). (3.21) 

On the other hand, the matter field contains only a conformally coupled mode and so we 
can use the exact propagator 

m^m^')) = ^^> (3-22) 



where Ax^ = x^ — x'^, Ax^ = r/^j^Ax^Aa;'". It should be noted that we have redefined the 
matter field. In fl2.20p and f l3.22p . the Schwinger-Keldysh indices are assigned as follows 

VAB = H^a{r)air')Axls, A, B = ±, (3.23) 

(3.24) 









Axl^ = 


-i\r- 


t' 


1 -ieY + {x 


-x')^ 








a4_ = 


-ir- 


1 

T 


+ ief + (x - 


x')^ 








Axl^ = 


-ir- 


r' 


- lef + (x - 


■x')^ 








Axl_ = 


-i\r- 


r' 


1 + ief + (x 


-x')^ 


where e 


is 


a positive 


infinitesimal quantity, 









From (I3.19P and (I3.2ip . the second term in (I3.18P is evaluated as 

\K'dMhn^{^){K)-^{^))dA{x)) ^^\oga{T){hl + y9^}0(x). (3.25) 

To evaluate the third term in (I3.18p . we need to perform the following integration 



-IK^d^j rfV CAB{^+{x)^A{x')){d,<P+{x)d'p<PB{x))d',<P{x') (3.26) 

iK^d^d^dpd„ I rfV CAB{'^+{x)LpA{x')){(t)+{x)(t)B{x'))ii){x'). 



Here we have not differentiated {(f{x)(f{x')) since it doesn't induce the IR logarithm. In 
order to evaluate the coefficient of log a(r), we may adopt the following approximation : 



t/V CAB{^+ix)VAix')){(t)+{x)(t)Bix'))(l){x') (3.27) 

{Mx)Mx)) [ d'x'[{Mx)Mx')) - {Mx)<l>-ix'))W). 



This approximation method has been introduced in Yukawa theory and scalar QED [T^ \TE\ 
m\. 

We argue that the IR logarithms appear only in the local terms. As we can observe in O2.20p . 
the logarithmic term in the propagator only depends on the relative coordinates. So they 
do not lead to an explicit time dependence unless in a singular limit, namely local terms. 
Therefore we extract local terms in quantum equation of motion. To do that, we need to 
expand up to the second order 

0(x') -^ 0(x) - d^^{x)Ax" + -do,dp^{x)Ax"Ax'^. (3.28) 

From ( K2T\f . ([322D, (ESID and dS^HD, (E^SD is written as 

^ log a(r)a,a.aA / ciV{ [^ - ^] 0(a:) - [^-^] 9.0(0:) (3.29) 

1 pAx^Ax^ Ax^Ax'^T^ o 7/ n1 



I 



■'++ 



To investigate the kinetic term, we need to perform the following integrations. We list the 
results below 

d<^d, j d^x' [^ - ^] = -M^5^H,\ (3.30) 

d,d,ds I d'x' [^ - ^] = 8^7rXW%', (3.31) 

d^dsdedr, J dx [-^^ A^2— J (3-32) 

We explain how to derive them in Appendix |Al In total, f l3.26p induces the following kinetic 
term with the IR logarithm 

^ log a{T){6^XdA + ^'^'d^ + ^\'d.d, + Sj'S^d.d^ + 6,\%d, + S^'d^d, 
+ KKW'd'o + S,\X5fd'}m. (3.33) 



From IKT9\i and IK33\i . the third term in fl318|) leads to 

^loga(r){-^9o^-^af}0(a:). (3.34) 

From (I3.25P and (I3.34p . the quantum equation of motion of a scalar field including the 
one-loop correction from soft gravitons is 

l + ^^logaiT))d^ix). (3.35) 

Although each contribution (I3.25p . (I3.34p breaks the Lorentz invariance, the total of them 
preserves it. The IR effect emerges just as an overall factor. Up to k^ if ^ log a(r), we can 
eliminate it by the following time dependent renormalization of a scalar field: 

^^Z^(P, Z^^l-^^loga(r). (3.36) 



We have checked that the same result is obtained in an exact calculation with the dimensional 
regularization. The IR logarithm originates from the dS symmetry breaking term in the 
graviton propagator in such a calculation. We also remark that the IR logarithm can be 
absorbed into the wave function renormalization factor even if we include a mass term as a 
perturbation. 

10 



Next we perform a parallel investigation with a Dirac field. The corresponding action is 

S= [ V^d^x tije^^Y^^ij, (3.37) 

where e'^g^ is a vierbein and 7*^ is the gamma matrix: 

7V + 7V = -2ri''K (3.38) 

The vierbein can be parametrized as 

e'^a = f^-V,, e^, = {e-^'r,. (3.39) 

3 

In a similar way to (I3.16p . we redefine the matter field Qiip — ). ip: 

S= I (fx itpe^^-f^Vf^ip. (3.40) 

By decomposing ip into the classical field and the quantum fluctuation 

^^^ + ^, (3.41) 

and differentiating (I3.40p with respect to ijj, the quantum equation of motion up to the 
one-loop level is written as 

tV^^rdMx) + t^{{h'^^UixW,)4x))rdMx) (3.42) 

Here we have approximated V^ ~ d^ since the differentiated gravitational field doesn't 
contribute to the IR logarithm. 

By substituting the identity {ip{x)%l){x')) = iriP^Ydp{4>{x)(f){x')) , 03.42p is written as 

tv'a^d.^i^) + ^^{{h^^)^{x){h^^)^{x))r^M^) (3.43) 



+ ^< / d'x' CAB{{h''a)4x){h\)A{x'))r7''{d,d,Mx)Mx')HdMx') ^ 0. 



From (I3.19P and (I3.2ip . the second term in (I3.43P is evaluated as 

k2 ^ k^H"^ 3 13 ^ 

^y((/^^p)4^)(/^''J-f(^))7"5.^(a:) ^ ^-^loga(r) x z{ - -7% + -j'd,}^{x] 

To evaluate the third term in (I3.43p . we need to perform the following integration 



K 



^ r]',J dV CAB{v+{x)Mx'))YY{dAMx)Mx')h'd,4^W) (3.45) 

'(¥.+ (x)^+(x)X7V7'5^5.5p I d'x'[{Mx)Mx')) - (0+(x)0_(x'))]V^(x'). 



11 



Here we have neglected the differentiated gravitational fields and adopted the same approx- 
imation procedure with scalar field case. 

Just like scalar field theory case, we need to evaluate local terms to estimate quantum IR 
effects due to soft gravitons. For such a purpose, we need to expand ^(x') up to the first 
order 

^{x') -^ ipix) - 9„^(x)Ax". (3.46) 

From ( Km . (1322]) and ( K^ . (EJS]) is written as 

'^^ \oga{T)r]P^ri'l%dudp f d^'x' (3.47) 



647r4 



VaxI 



1 - Ax" Ax" 



Ax^_-' '^ ^""^ '-Ax^^ Ax ^ 
By substituting ( 1330|) and ( 133T|) to ([MZD, (IMS]) is evaluated as 

^^loga(rX7'^7^7'{ - \iS,\% + d^d^d^ + Sj^Sfd,) + Is^W'do}^^). (3.48) 
From fl3:T9|) and (jMH]), the third term in KM\i is 

^loga(r) X ^{^l% - ^1%}^. (3.49) 

From (13.441) and (13.491) . the quantum equation of motion of a Dirac field including the one- 
loop correction from soft gravitons is 

(l + li^log«(^)) "" ^V^a^dJix). (3.50) 

Just like a scalar field, the IR effect from gravitons to a Dirac field preserves the Lorentz 
invariance. It can be eliminated by the following time dependent wave function renormal- 
ization of a Dirac field: 

7/;->Z^^, Z^ ^l-——^ log a(r). (3.51) 

We also remark again that the IR logarithm can be absorbed into the wave function renor- 
malization factor even if we include a mass term as a perturbation. 

We summarize our investigations in this subsection. Inside the cosmological horizon, the IR 
effects of the gravitons at the one-loop level preserve the Lorentz invariance both in a free 
scalar and Dirac field theory. We suspect that this is the case beyond the one-loop level in the 
both scalar and Dirac field cases. Although we think it is likely that soft graviton effects do 
not spoil Lorentz invariance of local physics, so far we have only demonstrated it by explicit 
calculations. We need to understand a mechanism to ensure it to all orders in perturbation 
theory. Our analysis has shown that the IR effects manifest as the overall factors of the 
kinetic terms. In free field theories, they can be renormalized away by a time dependent wave 
function renormalization. With interaction, the wave function renormalization contributes 
to the renormalization of the coupling constants. We investigate the IR effects in (f)'^ theory 
and Yukawa theory in the next subsection. 
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3.2 (f)^ theory and Yukawa theory 

Here we investigate the IR effects form gravitons to interacting field theories. As specific 
examples, we adopt 0^ theory and Yukawa theory. Since the coupling constants are dimen- 

3 

sionless, y/—g can be absorbed by the field redefinition Q(f) — >■ 0, Vt^ij^ ^ ip 

bU = -^0^ (3.52) 

6Cy = —g4>'4'4'- (3.53) 

After the wave function renormalization fl3.36p . f l3.5ip . the interaction terms are renormalized 

as 

5C, = -\K<\>\ (3.54) 

5£y = -gZ^Z^4njjil}. (3.55) 

First, we investigate 0^ theory. Up to the one- loop level, the following nonlinear terms should 
be added to the left-hand side in f l3.18p 

— Z'^^Hx) + ^ j d'x'd'x" CABCcD{{hn+{^){hnA{x')) (3.56) 

x(a^a,0+(x)0c(x"))(5X0B(x')0D(x"))0(a;')0'(x") 

+ ^ j d'x'd'x" CABCcDW)A{x'){hnc{x")) 

x(0+(a:)a;a>B(x'))(0+(x)a;a>^(x"))0(x)0(x')0(x"). 

Here we have performed the partial integration and neglected the differentiated gravitational 
field. The second and third terms denote the quantum correction to the vertex. The purpose 
in this subsection is to evaluate the effective coupling constant. To do that, we have only 
to extract the zeroth order of the classical fields in the Taylor expansion of the relative 
coordinates 

0(x')02(x"), 0(x)0(x )0(x") ^ 0^(x). (3.57) 

In a similar way to f l3.27p and fl3.45p . to evaluate these integrals up to ti^H^ log ^(t), we may 
adopt the following approximation 



((/i'^'^)+(x)(/i^'^) + (x)) / d^x'd^x". (3.58) 

In this approximation, only the following local term contributes to the remaining integrals 
(a^a,0+(x)0+(x')) = -(a^0+(x)9>+(a;')) ^ -i5^%^%^'\x - x'). (3.59) 
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As a result, f l3.56p is evaluated as 

- ^Z'^4>Hx) - \K'{ih''Ux){h''Ux))^'ix) (3.60) 

In the second line, we have substituted (13.191) . (13.211) and (I3.36p . (I3.60p implies that the 
effective coupling constant in 0^ theory decreases with cosmic expansion under the inffuence 
of soft gravitons 

Aeff ~ A{1 - — -^loga(r)}. (3.61) 

Next, we investigate Yukawa theory. Up to the one- loop level, the following nonlinear terms 
should be added to the left-hand side in (I3.42p 



-gZ^Z^<P{x)^{x) + ^— j (fx'dS" CABCcD{{h^a)+{x){h\)A{x')) (3.62) 

x(7'^9^7A+(x)V^C7(x"))(^D(x")9>B(x'))0(x")7'^(a;') 

+ i^ I d'x'd'x" CABCcD{{h''a) + {x){hnAix')) 

x{rd,Mx)Mx'')){d'^d'jB{x')Mx''mx'mx'') 

- t^jd'x'd'x" CABCcDma)A{x'){hnc{x")) 

x(^+(a;)9>B(x'))(0+(x)9;9>^(x"))0(a:")7>(x')- 

To evaluate the effective coupling constant, we have only to extract the zeroth order of the 
classical fields in the Taylor expansion of the relative coordinates 

0(x")^(x'), 0(x')^(x") -^ 0(x)^(x). (3.63) 

By substituting (^(x)-?/'(a;')) = iri'^^'j'^du{(f){x)(j){x')) and performing the parallel procedure 
with (13381) - (13:^9!) . (^m) is evaluated as 

- gZ^Zl^{x)i,ix) - ^((/,oo)4x)(/.°°) + (x))0(x)^(a:) (3.64) 

--^{1- i287r2 ^^Sa{T)}<Pix)i(){x). 

We find that the effective coupling constant decreases with cosmic expansion also in Yukawa 
theory: 

^cff ^ ^{1 -^^^^ log a(r)}. (3.65) 

As seen in (I3.6ip and (13.650 . the gravitational fluctuations outside the cosmological horizon 
influence the physics inside the cosmological horizon at the one-loop level. By the power 
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counting of the IR logarithms, the leading IR effect at the n-loop level is estimated as of order 
(K^iJ^ log a(r))". Since the perturbation theory is broken after k^ if ^ log a(r) ~ 1, we need 
a non-perturbative method to understand its long term consequences. Of course we need to 
make sure that this is a physical effect. In order to clarify such an issue, we investigate the 
gauge dependence of these perturbative IR effects in the next section. 



4 Gauge dependence 

In the previous sections, we have investigated the IR effects with the gauge fixing term fl2.12p . 
It is important to investigate the gauge dependence of the obtained results. In this section, 
we adopt the following gauge fixing term with a parameter /3: 

CoF = -^a^'-VF.F,, (4.1) 

F, = (3dph; - f3{D - 2)d^w + ^D - 2)h;dploga + ^{D - 2)wd,loga. 

This gauge fixing term coincides with (12.121) at /3 = 1. For simplicity, we consider the 
case /3^ — 1 ^ 1 since the deformation from fl2.12p can be investigated perturbatively. The 
deformation of the action at (9(/3^ — 1) is 

<5£^2_i ~ -i(/32 - l)a2 [ r7^-5^/i005,/i°° - 3doh""doh'' - -d.h^'^d.h'' (4.2) 

Here we have set -D = 4 and neglected massless conformally coupled modes. In addition, we 
have ignored ghost fields since they do not couple to matter fields. To investigate the gauge 
dependence, we evaluate the correction to the gravitational propagator from the additional 
term (g^D- 

As seen in (13.271) . (13.451) and (I3.58p . we only need to investigate the graviton propagator at 
the coincident point to evaluate the coefficient of k^ if ^ log a(r). From (I3.19P and (14. 2p . the 
IR logarithm can emerge in the following propagators 



{{h'')+{x){h'') + {x))\p._, ^ - z(/3^ - 1) y dV a\T')cAB (4.3) 

X { r/'-((/.o°) + (x)9;(/.°°)^(x'))((/.°°)+(a:)a:,(nB(x')) 
-3((/.°°)4^)5^(/^°°)A(a:0)((/^°°)+(^)5o(/^°°)B(^')) 

-H((^°°) + (^)^^(^°°)^(^'))((/^°°)+(^)5:(/^°°)i^(^'))}, 



{{h'')+{x){h^^)^{x))\p.., ^ - i{P' -I) j dV a\T')cAB (4.4) 

X { - ^((/^°°) + (x)9U/.°°)a(x'))((/^'^V(^)5;(/^'')b(x'))}, 
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{m^{x){h'^%{x))\p.., ^ - i{P^ - 1) y rfV a\T')cAB (4.5) 

x{((/^^^)4:^)5:n(/^'"^)A(x'))((/^'') + (:^)5;(/^'^^)B(x'))}. 



To evaluate them, we need to perforin the integration involving a massless and minimally 
coupled field 

-i{P'-l) f d'x' a'{T')cAB{Mx)dlMx')){M^)d',M^')). (4.6) 



In the integral, the following term of fl2.2Up contributes to the IR logarithm 

{ifixMx')) ~ --{ -logy + log a(r)a(r')} = --^logH^Ax''. (4.7) 

By substituting (14. 7p and using the following identity 

/-iVi;!^^-^^] =.-feMoga(r){«;C + i'w}. (4.8) 

the integral (14. 6 p is evaluated as 

-(/?' - l)^log«(r)«5,° + Ir/,.}. (4.9) 

We explain how to derive the identity (14. 8p in Appendix Rl 
From dM]), (i3])-(i3D are 

{ih'X{x){h'')4x))\,^-i ^ -iP' - 1) X -^^loga(r), (4.10) 



((n + (^)(/^^^)+(^))l/.-i=^0, 

H 
3' ' '^ 



{{K^)^{x)(h^')^{x))\f,._^ ^ -(/32 - 1) X {^H^' + ^'b^"" - -5^^5^')_ log a(r). 



We should note that the deformation of propagator is proportional to the original one. As a 
result, in the deformed gauge ( 14. ip . we have only to replace the gravitational propagator as 
follows 

{{hn+{x){hn+{x)) ^ {1 - (/?^ - i)]{{hn+{x){hn+{x)). (4.ii) 

From this fact, we can conclude that the Lorentz invar iance is preserved for a continuous /3 

Z,^l-{2-P'f-^\oga{T\ (4.12) 

Z^^l-(2-/52)|^loga(r). (4.13) 
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However the effective couplings are found to be gauge dependent. 



Aeff ^ A{1 - (2 - /3')^^0^ loga(r)}, (4.14) 

^eff-^{l-(2-/3^)^^loga(r)}. 



Here we summarize our findings as follows. In the original gauge the dimensionless coupling 
constants A and g are screened by soft gravitons and decrease with time in de Sitter space. 
They acquire nontrivial scaling exponents in de Sitter space as 

21 '^Q 

Aeff = A/(t)X, g^s = gf(t)ir. (4.15) 

where f{t) = 1 — ^-^rHt. In order to confirm that this is a physical effect, we have examined 
the gauge parameter dependence of these results using one parameter family of the graviton 
propagator. As it turns out f{t) depends on a gauge parameter /3 as f{t) = l — {2—j3'^)^^j^Ht. 
Neverthless we observe that the following relative scaling relation is independent of a gauge 
parameter. 



^cff~(Aeff)56 (4.16) 



We point out that the situation here is analogous to the scaling exponents of the operators 
in two dimensional quantum gravity. The scaling exponents of the individual operators are 
gauge dependent. However the relative scaling exponents are gauge invariant [I8],|19]. It is 
because there is no unique way to specify the scale there. Analogously in our case there is no 
unique way to specify the time as it depends on an observer. Here a sensible strategy is to 
pick a particular coupling and use its time evolution as a physical time. In this setting the 
relative scaling exponents measure the speed of the time evolution of the couplings in terms 
of a physical time. Although the choice of time is not unique, the relative scaling exponents 
are gauge independent and well defined. We still need to check the validity of this picture 
against large deformations of a gauge parameter. 



5 Conclusion 



The inflation theory postulates that almost scale invariant density perturbation is generated 
by quantum fluctuations in a de Sitter type space. The detection of almost scale invariant 
gravitational fluctuations is a crucial test of the inflation theory since it universally generates 
them. Thus the propagator of gravitons contains de Sitter symmetry breaking IR logarithms 
in de Sitter space. 

In this paper we have investigated its physical implications on microscopic physics. Namely 
we have investigated the effects of super-horizon modes of gravitons on the dynamics of the 
sub-horizon modes of matter fields. By evaluating the kinetic terms of scalar and Dirac fields 
up to the one-loop level, we have found that the IR effects from gravitons preserve the Lorentz 
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invariance. In particular the velocity of massless particles remains universal irrespective of 
the spins. Furthermore they can be absorbed into the wave function renormalization of 
scalar and Dirac fields. S. B. Giddings and M. S. Sloth investigate the scalar two point 
function in the quenched approximation [20]. S. P. Miao and R. P. Woodard show that the 
IR effect on the Dirac field is recognized as the wave function renormalization |13j. Their 
results qualitatively agree with ours. 

In the interacting field theory with dimensionless couplings, we find that the couplings be- 
come time dependent and scale with time. In such a situation any coupling can be chosen 
as a physical time. Then the relative scaling exponents of the couplings measure the speed 
of time variation of the couplings in terms of a physical time. In fact we have checked that 
they are invariant under an infinitesimal change of gauge. 

The results obtained in this paper are the one-loop effects. Since the IR effects at each loop 
level manifest as polynomials in log a(r), perturbation theories are broken after enough 
time passed. In other words, the obtained results describe physics at the initial stage 
K^ if ^ log a(r) ^ 1. To investigate the eventual contributions to physical quantities, we 
need to evaluate the IR effects nonperturbatively. Although the IR effects from specific mat- 
ter fields have been investigated nonperturbatively [U [TOl [161 [13 [H] ) the nonperturbative 
approach for the IR effects from gravitons is an open issue. 

Since our Universe is dominated by dark energy, our results imply the couplings of the 
standard model become time dependent. Since the coefficients of IR logarithms are of 
0{k'^H'^) ~ 10~^^°, it is not observable now. However it is relevant to the ultimate fate 
of the Universe as their effect grows linearly with cosmic time like Hawking radiation. This 
effect is much larger in the inflationary era as k^H^ could be as large as 10~^. 

Since the cosmological constant is a function of the couplings of the microscopic theory, it 
may acquire time dependence if the couplings evolve with time. So we need to investigate 
such effects to understand possible time dependence of the cosmological constant in addi- 
tion to pure matter and gravity contributions separately. Let us follow the change of the 
couplings and cosmological constant under time evolution. We note here that the vanishing 
cosmological constant is a fixed point of the evolution where the couplings become constant. 
It seems conceivable that such a self-tuning mechanism may play an important role in the 
cosmological constant problem. 
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A Derivation of ( 13:301 ) -( 13:321 ) and (\4M 



Here we explain how to derive the identities ( I3.30p -f l3.32p and (14. Sp . Each integrand in their 
left-hand sides is written as 



(A.i: 



1 

Ax2 


= ^dHogH^Ax\ 


Ax„ 

A „9 


= ldo.\ogH^Ax\ 



(A.2) 



^^^ = ld^{Ax,\ogH'Ax') - Iva.logH'Ax', 



(A.3) 



AXaAX/3 1 



Ax^ 4;-"-M 2 



{d^d0--Vapd'}\ogH'A 



X 



(A.4) 



where we abbreviate the indexes ++, H — because the above identities work in both cases. 
By using them and leaving differential operators out of the integrals, the left-hand sides in 
(E30])-(E32D and gSD are 



= - ^5,V^o /^Vflogif^Ax^^ - \ogH'Axl_], 



(A.5) 



dfjd^ds 
1 



I d\'[ 



AXr 



AXr 



Ax\^ Ax\. 



-] 



K'^pW'dt / d'x'[\ogH'Axl^-\ogH'Axl_\ 



(A.6) 



a a n n /" J4 ,[^^aAxp AXaAxii. 

^dsdsdsjdx[^^-^^] 

- - l^^WWXd'o I d'x' Arflogff^Ax^^ - logi:f^Ax^_] 
- ^v^p6^%\Xd'o J d'x' [ log H'Axl^ - log H^Axl_] , 



(A.7) 



4 , 1 rAxaAx/3 Aa;aAx/3. 

^^ 7;^[^:;4 ^:::^] 



r'2L Ax^ 



Axi 



^{CV + \v^pH I d'x' ^^[XogH^Axl^ - logiJ^Ax^_] , 



(A.8) 
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Note that the differential operator outside the integral is equal to the time derivative 8^ — ?■ 
S^Oq. Furthermore, we have replaced Ax^ to —5^/S.t in (]A.7p . 



As a concrete example, we calculate the following integral 



dl I d^x' [ log H^Axl^ - log H^Axl_] . 



From f l3.24p . the logarithm with each index ++, H — is 



(A.9) 



-2 2\ 



logH^Axl^ = logH'lAr^ - r'\ + inOiAr' - r 

log/f^Ax^. = logif^lAr^ - r^\ - ine{AT^ - r^){e{Ar) - 0{-At)}, 

where r^ = (x - x')^. By substituting ( lA.lOl) . ( JA.91) is 



(A.IO) 



d^ / dr' / Anr^dr 2m 



nn 



d^ / dr' At^ = leivrl 
3 ./_ 1 



In a similar way, the other integrals are evaluated as 



^0^ /"dV AT[\ogH^Axl+ - \ogH^Axl_] = Um^ 



(A.ll) 



(A.12) 



dl I d^x' \ [ log if 2Ax^+ - log H^Axl_] = 16z7r2 T dr' % (A.13) 



IGivr^ log a(r). 



We have extracted the IR logarithm in the last line of ( 1A.13I) . By substituting ( ]A.llj) - (lA.13l) 
to flA.SP - flA.Sp . we obtain the desired identities 



dc^dp j d^x'[—, 



2 A„2 



-] = -4.vr^CV> 



(A.14) 



dpd-yds I d X 



j4 J r ^^a /AXc 



LAx^+ Ax^. 



-] = 8^7rX%,%^%s', 



(A.15) 



d^/dsdedr, / (i'^x'f 



Axl_ 



Ax^ 



Ox 0;r Ox 



- 32z7r^5,"<5;5,"<5,"<5/V - 8m^r^.,5;5,"5,'^V' 



(A.16) 



^^-' i[^^ - ^S:^] - -4-^loga(r){VC + \n.A- 



r- 



Axt 



Axi 



(A.17) 
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